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Estimates of the bulk metal abundance of the Sun derived from the latest 
generation of model atmospheres are significantly lower than the earlier standard 
values. In Paper I we demonstrated that a low solar metallicity is inconsistent 
with helioseismology if the quoted errors in the atmospheres models (of order 
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. 0.05 dex) are correct. In this paper we undertake a critical analysis of the solar 

^ . metallicity and its uncertainty from a model atmospheres perspective, focusing 

O , on CNO. We argue that the non-LTE corrections for abundances derived from 

O ■ atomic features are overestimated in the recent abundance studies, while sys- 

tematic errors in the absolute abundances are underestimated. If we adopt the 
internal consistency between different indicators as a measure of goodness of fit, 
we obtain intermediate abundances [C/H] = 8.44 +/- 0.06, [N/H] = 7.96 +/- 0.10 
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c/2 ! and [0/H] = 8.75 +/- 0.08. The errors are too large to conclude that there is a 

^ ' solar abundance problem, and permit both the high and low scales. However, the 



X 



center-to-limb continuum flux variations predicted in the simulations appear to 
be inconsistent with solar data, which would favor the traditional thermal struc- 
', ture and lead to high CNO abundances of (8.52, 7.96, 8.80) close to the seismic 

scale. We argue that further empirical tests of non-LTE corrections and the ther- 
mal structure are required for precise absolute abundances. The implications for 
beryllium depletion and possible sources of error in the numerical simulations are 
discussed. 

Subject headings: convection, Sun:abundances, Sun:evolution, stars: abundances 



1. Introduction 



The uncertainty in the absolute chemical composition of stars is the limiting factor in 
our ability to do high precision stellar astrophysics. Traditionally, we have had to rely on a 



- 2 - 



small database of fundamental stellar parameters such as mass, distance, and radius. How- 
ever, current and upcoming space missions promise a wealth of astromctric and photometric 
data. Large surveys undertaken primarily for other purposes (microlensing, planet searches 
and cosmology) have discovered thousands of eclipsing binaries, yielding numerous precise 
mass estimates. The rapidly developing field of optical interferometry has also permitted a 
growing number of direct radius estimates. Asteroseismology is also growing in importance, 
and missions such as COROT promise a wealth of detailed information on the pulsational 
properties of solar-like stars. 

Our stellar interiors models have become highly sophisticated and successful when com- 
pared with observational diagnostics. In particular, the resolution of the solar neutrino prob- 
lem in favor of the solar model predictions and the agreement between theoretical predictions 
and helioseismic data are both encouraging signs. The combination of better observations 
and theory has opened the prospect of a new era of precision stellar astrophysics, which 
could have broad consequences for diverse subfields of astronomy. 

Stellar atmospheres theory has traditionally employed a series of approximations when 
deriving abundances. Classical models assume an ad hoc turbulent velocity field adjusted to 
yield abundances independent of excitation potential and line strength. Convection is usually 

treated in an approximate fashion, with the mixing length theory. Horizontal temperature 
fiuctuations (granulation) are not inchidcd. The models also typically assume that the 
molecular and atomic levels arc described by local thermodynamic equilibrium (LTE), e.g. 
by the local temperature alone. The compilations of solar abundances used for theoretical 
solar models (Anders & Grevesse 1989; Grevesse & Noels 1993; Grevesse & Sauval 1998) 
employed model atmospheres with approximations at the level described above. The mean 
thermal structure employed was semi-empirical (Holweger & Mueller 1974, hereafter HM74). 

When these approximations are relaxed, different conclusions about the abundances are 
obtained. Departures from LTE are expected at a modest level for solar conditions, and have 
been investigated by a number of authors (for example Carlsson 1986; Stuerenburg & Hol- 
weger 1990; Kiselman 1993; Shchukina & Trujillo Bueno 2001; Wedemeyer 2001). Numerical 
simulations of convection have matured to the level where they can be used to predict veloc- 
ity fields, temperature fluctuations, and changes in the mean thermal structure of the upper 
atmosphere (Stein & Nordlund 1998). Abundances derived from these simulations yield a 
very different pattern, which has been developed in a series of papers (Allendc Prieto, Lam- 
bert & Asplund 2001, 2002; Asplund et al. 2000a,b, 2004, 2005); papers by Lodders (2003) 
and Asplund, Grevesse & Sauval (2005) summarize the revised abundance scale. 

The net effect is in the sense of systematically lower metal abundances. The downward 
revisions for the heavier elements (e.g. Fe and Si) are small, while the claimed reduction 
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in the abundances of lighter species (especially CNO) is more dramatic. Models employing 
different treatments of granulation and non-LTE corrections (Holweger 2001; Steffen & Hol- 
weger 2002) predict smaller abundance reductions. The central temperature predicted by 
interiors models is sensitive to the abundances of the heavier elements, but not the lighter 
ones. As a result, the new abundance scale docs not disturb the agreement between inte- 
riors models and observational data for purposes such as the mass-luminosity relationship 
and solar neutrino fluxes. However, the inferred solar sound speed profile, and the radii of 
interiors models, is sensitive to the bulk metallicity. Serious problems have emerged when 
comparing interiors models with the revised abundance scale. These discrepancies are ev- 
idence for problems in our understanding of stellar interiors, stellar atmospheres, or both. 
In Paper I (Delahayc & Pinsonneault 2006) we investigated the errors in solar abundances 
predicted by the combination of stellar interiors models and helioseismic data. In this paper 
we examine the uncertainties in the abundance predictions from stellar atmospheres theory. 
Wc begin with a brief summary of the results from Paper 1, and follow with a discussion of 
the motivation and main results from the revised stellar atmospheres models. In Section 2, 
we perform a critical analysis of the precision of the solar CNO abundances and discuss the 
implications for Be. We demonstrate in that section that the errors in the abundances are 
larger than previously estimated, and that there is evidence that the "best" current solar 
CNO abundances are intermediate between the new and old scales, with errors permitting 
both. We discuss the implications of our finding and future tests in Section 3. In particular, 
we argue that inconsistencies between the solar thermal structure and that predicted by the 
simulations would favor a higher abundance scale closer to the seismic value and discuss 
uncertainties in the numerical convection simulations. 



1.1. Constraints from Helioseismology 

Helioseismology provides two powerful constraints on the solar composition: diagnostics 
of the internal solar temperature gradient and diagnostics of the equation of state. Inversions 
of the observed solar pulsation frequencies yield accurate measures of the sound speed as 
a function of depth. In turn, the gradient in the sound speed can be directly tied to the 
temperature gradient. Since the temperature gradient is related to the opacity, and thus 
the composition, information on the solar abundances is encoded in the seismic data for the 
radiative interior. One can even obtain meaningful constraints on the solar age from the 
helium abundance profile deduced in the deep interior. 

In addition to the vector information on the sound speed profile, there are also precise 
scalar quantities that can be extracted. The thermal structure at the base of the solar 
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convection zone is nearly adiabatic, while the temperature gradient in the interior is radiative. 
As discussed in Paper I, Section 2.2 the resulting discontinuity in V generates a distinct signal 
that can be used to precisely localize the base of the solar convection zone (i?c2 = 0.7133 ± 
0.0005Rsun Basu & Antia 2004). Seismology also sets strict limits on convective overshooting 
(< 0.05Hp, Christensen-Dalsgaard et al. 1995). The depth of the solar convection zone is 
sensitive to the light metal abundances in the Sun but insensitive to most of the other 
uncertainties in solar interiors models (see Paper I for a detailed error budget). Ionization 
also induces a depression in the adiabatic temperature gradient, and the absolute abundances 
of the species in question can be inferred from the magnitude of the perturbation in the 
surface convection zone. An extremely precise surface helium abundance can be deduced 
from this effect {Ygurf = 0.2483 ± 0.0046, see Paper I, Section 2.3 for the sources used in 
this estimate). More recently, Antia & Basu (2006) have demonstrated that the ionization 
signal of metals in the convection zone can be detected in the seismic data, leading to a 
bulk metallicity Z = 0.017 ± 0.002. Because the majority of the solar metals are in the 
form of CNO, this is primarily a constraint on their abundance. In principle one might be 
able to use this technique to solve for individual heavy element abundances by fitting the 
strength of distinct ionization stages. However, it is not yet clear that there is sufficient 
spatial resolution in the seismic data to permit such a detailed analysis. 

In Paper I we demonstrated that the combination of the surface convection zone depth 
and surface helium abundance constraints was a powerful diagnostic of the solar heavy ele- 
ment abundances. The surface helium abundance is tied to the initial solar helium abundance 
with a correction for gravitational settling. The initial helium is sensitive to the central opac- 
ity and the abundances of the heavier metals (especially iron). The convection zone depth is 
sensitive to the opacity at temperatures 2 million K, where bound-free opacity from light 
metals (CNONe) is an important contributor. The most significant new finding in Paper I 
was that the combination of the two scalar constraints could be used to rule out some abun- 
dance combinations with high statistical significance. The detailed sound speed profile adds 
additional information; we found that models consistent with the scalar constraints could 
be constructed with low oxygen and very high neon, but such models exhibited substantial 
sound speed deviations relative to solar data in the deep interior. 

The inferred solar oxygen and iron abundances {[0/H] ~ 8.86 ± 0.05, [Fe/H] — 7.50 ± 
0.05) are consistent with the Grevesse & Sauval (1998) absolute abundances, but strongly 
inconsistent with the new abundance scale within its quoted errors (Lodders 2003; Asplund, 
Grevesse & Sauval 2005). Although there are potentially positive chemical evolution conse- 
quences for the revised abundance scale (Turck-Chieze et al. 2004), it is not easy to generate 
interiors models that are consistent with both seismology and the low abundance scale. The 
most commonly cited possible explanations on the interiors side (high neon, enhanced grav- 
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itational settling, and errors in the high temperature radiative opacities) are all strongly 
disfavored. As previously mentioned, solutions with high neon degrade agreement with the 
sound speed profile, and are also problematic from a stellar atmospheres perspective (Schmelz 
et al. 2005). An increase in the degree of gravitational settling increases the convection zone 
depth but decreases surface helium, trading improved agreement with one diagnostic for 
worse agreement in another. Enhanced differential settling of metals with respect to helium 
is inconsistent with the underlying physics and would have to be extreme (Guzik et al. 2005). 
Three independent quantum mechanical calculations yield extremely similar Rosseland mean 
opacities at the temperatures of interest for the base of the solar convection zone. As dis- 
cussed in Paper I, both the atomic physics and equation of state are relatively simple in this 
regime, and the concordance between different calculations is thus not a surprise. Physical 
processes neglected in classical stellar models (such as rotational mixing and radiative accel- 
eration) can be independently constrained by other data, and in any case they would tend to 
induce higher rather than lower surface abundances. The scalar constraints are insensitive 
to the other theoretical ingredients in standard solar models (e.g. convection theory, sur- 
face boundary conditions, low temperature opacities, equation of state, and nuclear reaction 
rates) . The considerations above indicate that it is extremely challenging to reconcile a low 
solar metal abundance with current stellar interiors models and seismic data. This does not 
imply that a metal-poor Sun is impossible, but it certainly motivates an investigation of the 
uncertainty in the atmospheres models used to derive the abundances. 



1.2. Model Atmosphere Ingredients 

The new solar abundance estimates are derived from a variety of changes in the model 
atmospheres. Changes in oscillator strengths and equivalent widths of spectral lines con- 
tribute for some diagnostics, and as discussed below we largely concur with the revised 
values. The magnitude of non-LTE corrections depends on the atomic model and the rel- 
ative importance of photo-excitation and collisions on the level populations in the model 
atmosphere. The comprehensive re-examinations of the solar oxygen (Asplund et al. 2004, 
hereafter AGSAK04) and carbon (Asplund et al. 2005, hereafter AGSAB05) adopted a partic- 
ular model for NLTE corrections, and we assess its accuracy and uncertainties by comparison 
with limb darkening data and other published calculations. In AUende Prieto, Lambert & 
Asplund (2001, 2002) the abundances derived from forbidden lines have been reduced by the 
application of blending corrections; AGSAK04 and AGSAB05 used the revised equivalent 
widths for C and O abundance studies. Both the uncertainty in these corrections and their 
central value is incorporated in our error analysis. 
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Three other coupled changes in the atmosphere models are the treatment of convec- 
tive velocity fields ("macro/microturbulence"), horizontal temperature fluctuations (granu- 
lation), and the impact of convcctive overshooting on the mean thermal stratification. All 
of these features are derived from numerical convection simulations; a good discussion can 
be found in Stein & Nordlund (1998). Their combined impact can be deduced by the com- 
parison of the results in the 3D case in the published studies with the results from the 
semi-empirical ID HM74 thermal structure and the theoretical ID MARCS models. Since 
the initial 3D model atmosphere is derived with physics similar to that in the MARCS code, 
the impact of the convection treatment can be indirectly inferred by comparing MARCS and 
3D abundances. A comparison of HM74 with MARCS and 3D is a measure of the impact 
of different choices of the thermal structure. 

The numerical convection simulations predict line profiles in excellent agreement with 
the data for iron and sihcon hnes (Asplund et al. 2000a,b). There are some trends with 
excitation potential that may be related to NLTE corrections (Shchukina & Trujillo Bueno 
2001) and issues with the initial generation of simulations when compared with line profiles 
in the outer solar photosphere (Scott et al. 2006). The concordance of the predicted ampli- 
tude of horizontal temperature fluctuations with the solar granulation pattern is encouraging 
(Stein & Nordlund 1998; Asplund et al. 2000a). The validity of the steep solar tempera- 
ture gradient in the simulations is not as clear; there is an apparent conflict between the 
simulations and the mean thermal structure of the solar atmosphere, as well as the degree 
of convective penetration for the upper atmospheric layers (Ayres et al. 2006). Since these 
effects are all tied together in the abundance studies, we focus on the agreement between 
different diagnostics of abundance as a valuable test of the precision of the results obtained 
from different atmospheric models. 



2. The Solar CNOBe Abundances 

Revised solar abundances have been derived for a number of species. In this paper 
we focus on CNO for several reasons. First, these are the elements where the difference in 

abundance is largest, and they are also the cases where there is the largest variety of distinct 
abundance indicators. The difference between interiors and atmospheres based abundances 
for heavier species, such as Fe and Si, are not statistically significant. Furthermore, details 
of the new abundance estimates and internally consistent comparison with prior work are 
published for only some of the heavier elements. The solar Be abundance is an important 
diagnostic of mixing, and the photospheric abundance is linked to the solar O (Balachandran 
& Bell 1998). We therefore also briefly discuss the imphcations of our result for Be. We 
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begin with a description of our overall approach, and then follow with individual sections 
on oxygen, carbon, and nitrogen. Our primary references for O, C, N are respectively 
AGSAK04, AGSAB05, and Asplund, Grevesse & Sauval (2005). We include other studies 
for external comparisons, and discuss newer results based on other models and diagnostics 
when available. 



2.1. Overall Approach 

Our primary metric for the accuracy of the abundances derived from the different model 
atmospheres is the consistency of the estimates derived from distinct classes of indicators. 
One might justifiably apply a different standard, noting for example the greater degree of 
sophistication in the input physics for the hydrodynamic simulations. However, it is not clear 
that the ingredients that induce the abundance changes (such as changes in the temperature 
gradient) are actually a necessary consequence of these improvements in the model, as the 
absolute errors in the first-principles theoretical models have not been quantified. We proceed 
in two steps. In the first, we compute relative abundances and errors within a given assumed 
model atmosphere for the atomic and molecular features. Since the abundances from atomic 
lines have smaller model-to-model differences (and thus smaller systematic errors), we adopt 
the atomic abundances for our base estimate. The difference between atomic and molecular 
abundances is then used to infer which of the different atomic scales should be adopted for 
the central value, and the uncertainty in the differential scales is used as a measure of the 
systematic error arising from the choice of model atmospheres. 

For the atomic line diagnostics, we include random errors from the dispersion of results 
from single permitted lines about the adopted mean. For the forbidden lines, uncertainties 
in oscillator strengths and blending corrections become the dominant random error source. 
We also include systematic errors (NLTE corrections and zero-point shifts in the average 
oscillator strength) by comparing the study values with external constraints and other pub- 
lished calculations. This is important when comparing the atomic and molecular indicators, 
since NLTE corrections are usually included in the former but not the latter. As a result, 
changes in the degree of NLTE corrections have a direct impact on goodness of fit. When 
available, we adopt a weighted mean of the permitted and forbidden atomic indicators and 
the error in the mean when comparing with the molecular data. 

For the molecular indicators, we include all measured lines of a given diagnostic and use 
the total dispersion (rather than the error in the mean) as a measure of goodness of fit. We 
adopt a weighted mean of the various indicators for an average molecular abundance. How- 
ever, the error in the mean will understate the true uncertainty; it is frequently the case that 
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the mean values for different molecular diagnostics differ by more than the dispersion within 
each indicator. We therefore compute the dispersion of each indicator about the adopted 
mean and average these values to obtain an error in the molecular abundances. One could 
alternately compute the dispersion in the molecular abundances and treat this as a measure 
of the systematic uncertainty, adding it to the error in the mean in quadrature; this proce- 
dure yields somewhat smaller errors. Although the latter approach may be practical when 
there are numerous molecular probes available, we prefer the former method for situations 
like oxygen (where there are 2 values, and thus an unreliable estimate of the uncertainty in 
the mean). 

We derive final abundance estimates for each species by comparing the mean abundances 
derived within each class of models for atomic and molecular indicators. In the case of 
oxygen, the 3D and HM74 models exhibit comparable differences with opposite sign, while 
the MARCS models have an internally consistent intermediate abundance. We therefore 
adopt a mean of the different derived oxygen abundances and an uncertainty from the scatter. 
In the case of nitrogen, the HM74 abundances are preferred. The case of carbon depends on 
the origin of systematic differences in abundances inferred from atomic features. If the low 
zero point of AGSAB05 is adopted, the 3D models are the favored solution. If the higher 
zero point of previous work is adopted, the situation is similar to that for oxygen. 



2.2. Oxygen Abundance Indicators 

AGSAK04 derived a low solar oxygen abundance (8.66 ±0.05) from four distinct indica- 
tors: atomic lines (forbidden and permitted) and two different classes of infrared molecular 
lines ((v,r) and (r,r)). All four indicators had formerly been used to obtain higher abso- 
lute oxygen abundance (8.83 to 8.90). Abundance estimates from all indicators are reduced 
in the theoretical atmospheres that include substantial overshooting because fines become 
stronger for a fixed abundance in the presence of a steeper temperature gradient. Molecular 
abundances are reduced more than atomic ones because the cooler atmospheric structure of 
the 3D hydro models changes the chemical equilibrium. AGSAK04 included other effects 
that reduced the abundances derived from atomic features without impacting the molecu- 
lar indicators: a combination of changes in oscillator strengths, the inclusion of blending 
features, and their claim of large non-LTE effects for the permitted lines. 

In this section we discuss the uncertainties in each of these cases. We advocate a 
substantial decrease in the magnitude of the NLTE corrections for the [0/H] derived from 
the permitted atomic oxygen lines, and an increased error in both the abundance derived 
from the forbidden line (from uncertainties in the oscillator strength) and the permitted 
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lines (from uncertainties in the NLTE corrections). We also derive an increased error from 
the [0/H] derived from the IR molecular lines from the internal scatter and trends with 
excitation potential, and argue that the correspondence between the trends in the MARCS 
and 3D models is evidence for errors in the common underlying model. Our basic results 
are summarized in Table 1. In Table 1, the upper part of the table repeats the mean 
abundances and errors for the 4 different indicators presented in AGSAK04. We present our 
revised estimates for the same three cases in the lower part. The last three rows in each 
sub-table give the mean discrepancy between atomic and molecular indicators. At the end 
of the section we synthesize this information to obtain our best estimate for the solar oxygen 
abundance. In each of the following subsections, the error estimates derived are internal 
ones. The differences between the results from the three classes of models are prima facia 
evidence that systematic errors are important. The systematic errors are discussed in the 
final subsection. 



2.2.1. Forbidden Oxygen Lines 

The forbidden oxygen line at 6300.3 A has traditionally yielded high oxygen abundances. 
AUende Prieto, Lambert & Asplund (2001) argued that a nearby blended Ni line contributed 
significantly to the oxygen feature. They treated the continuum level, log (gfNi), and the 
oxygen abundance as free parameters, but assumed that the line profiles as given from the 
simulations were exact. The inclusion of the Ni feature induced a direct reduction of 0.13 
dex in the inferred oxygen abundance. In addition, the usage of a 3D model atmosphere 
structure led to a further reduction of 0.08 dex in the oxygen abundance to an estimated 
[0/H] = 8.69 +/- 0.05. In AGSAK04 another forbidden hue at 6363.7 A was considered as 
a second indicator. The authors reduced the equivalent width by 0.5 mA for an estimated 
contribution from a blended CN feature to obtain an oxygen equivalent width of 1.4 mA , also 
implying a low abundance. We begin our analysis by noting that abundances derived from 
blended features are usually treated with great caution. The most conservative procedure 
is to ignore the blending feature and treat the derived abundance as an upper limit. When 
this is done for the two forbidden lines, the maximum abundance obtained for 3D, HM, and 
MARCS are (8.82, 8.8), (8.9, 8.88), (8.86, 8.84) respectively. In the section that follows, we 
include the reduction in abundances from estimates of the blending contribution. 

We adopt the AGSAK04 values for [0/H] derived from the 6300.3 A line, subject to the 
caution on the strength of the Ni blending feature below. For the 6363.7 A line, Melcndez 
(2004) argued that the (10,5) Q2 25.5 CN line is unblended in the solar spectrum and has 
the same oscillator strength as the feature blended with the forbidden hue. He derived a 
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smaller correction for the blended CN line (0.35 mA rather than the 0.5 mA value used in 
AGSAK04), which we adopt here. This leads to a modest 0.04 dex increase in the derived 
abundance from that line, which we also treat as an uncertainty in the [0/H] derived from 
this feature from the uncertainty in the contribution of CN to the blended feature. The 
error analysis for a blended feature is more complex than the one that can be employed 
for an isolated line. The derived [0/H] is sensitive to the continuum level, and an error 
component for this should be included; AUende Prieto, Lambert & Asplund (2001) estimate 
this uncertainty at 0.02 dex, which we adopt for both lines. AGSAK04 adopted lower 
values for the oscillator strengths than those found in the NIST database, but their choice 
is well-supported by the improved atomic physics (see Storey Sz Zeippen 2000) . However, 
the errors in individual theoretical log(gf) values are higher than those assigned in AUende 
Prieto, Lambert & Asplund (2001) and subsequent papers; we adopt 0.04 dex for individual 
lines. AUende Prieto, Lambert & Asplund (2001) also estimated that uncertainties in the 
underlying equation of state induces errors of 0.02 dex. 

Especially for the 6300 A line, the results depend heavily on the detailed line profiles, 
particularly in cases where the individual components cannot be directly disentangled. Al- 
lende Prieto, Lambert & Asplund (2001) estimated uncertainties of 0.02 dex from the central 
wavelength of the Ni feature and a 0.04 dex uncertainty from the central wavelength of the 
[O I] hne. We treat these errors as representative of the uncertainties in the hue profiles, 
and adopt them for both forbidden lines. 

The treatment of the Ni line in the main forbidden line is more problematic. In the initial 
study, the oscillator strength was highly uncertain. AUende Prieto, Lambert & Asplund 
(2001) treated log (gf) for Ni as a free parameter. The continuum level, log{gfMi)-, and 
[0/H] were treated as free parameters and the combination that produced the minimum 
was adopted. However, Johansson et al. (2003) have measured the oscillator strength of the 
Ni feature (log gf = -2.11), and the Ni abundance of the solar mixture is well-constrained by 
knowledge of the solar Si/Fe and the relative mcteoritic abundances. With the new gf value 
and [Ni/H]=6.25 the blending feature would be 0.23 dex stronger than the best value 
obtained in the 2001 paper. AGSAK04 did not report the inferred strength of the blending 
feature in their fit, but the similarity in the absolute abundance suggests that it would be 
comparable. We believe that it is no longer appropriate to treat this as a free parameter, 
and the same method should be used as is done for other blended features: namely, the 
strength of the Ni line should be held fixed (and varied within its uncertainty) while the free 
parameters are the oxygen abundance and continuum level. 

The direct effect of increasing a blending contribution is usually to decrease the abun- 
dance, but the present case is more comphcated. For the quoted values in the AGSAK04 
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fit, the Ni line contributes 25% of the total equivalent width of the line. An increase of 0.23 
dex in the strength of the feature would imply a total Ni contribution of 43% of the blended 
equivalent width at a fixed continuum level; such a combination would be a poor fit to the 
line shape and would yield a reduction in the inferred oxygen of 0.12 dex. An increase in 
the continuum level would be required to restore the agreement with the line profile, which 
would in turn lead to an increased total equivalent width. The net impact on the derived 
oxygen is not obvious, and not necessarily in the negative direction. Reetz (1998) invoked 
log gf=-1.95 for the Ni feature and obtained [0/H]=8.75 for the forbidden hne when his 
oscillator strength for the forbidden hne was adjusted to the same value as that employed 
by AUende Prieto et al. In the absence of other information, we assign an additional error 
component of 0.04 dex for the strength of the Ni feature, slightly higher than the value 
advocated by Melendez (2004). 

Adding the errors in quadrature, we obtain an uncertainty of 0.078 per line (0.055 in 
the average) and abundances derived from the forbidden lines systematically 0.02 dex higher 
than those found in AGSAK04. 



2.2.2. Permitted Atomic Oxygen Lines 

The permitted atomic oxygen lines have relatively high oscillator strengths, but very 
high excitation potentials. The most commonly used features are the 01 triplet at 7771.8, 
7774.2, 7775.4 A , with an excitation potential of 9.15 cV. AGSAK04 also considered three 
other atomic features (6158.1 A , 8446.7 A , and 9266 A ). The primary reason for the low 
oxygen abundance inferred by AGSAK04 from the triplet is a large non-LTE correction. 
Because these lines arise from such a high energy state, non-LTE effects must be included. 
However, the quoted values of the non-LTE corrections in the literature vary drastically. 
For the triplet, the average is -0.06 dex for Holweger (2001), -0.22 to -0.28 for AGSAK04, 
and -0.16 dex for AUende Prieto et al. (2004). These variations can be partially traced to 
different assumptions about the importance of coUisional excitation (as opposed to photoion- 
ization), but larger differences at the 0.10 dex level remain even for cases that make similar 
assumptions about hydrogen collisions. 

AGSAK04 neglected hydrogen collisions in their estimate of the non-LTE effects. They 
justified this by noting that for some well-studied hues, the classical Drawin (1968) formuhsm 
overestimates the collision rate. However, an inspection of their Figure 6 indicates that the 
neglect of collisional effects in their model yields changes as a function of limb darkening 
that differ from the observed solar values. This impression is confirmed by the more detailed 
study of AUende Prieto et al. (2004), who found that models including hydrogen collisions 
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( their Sh — I case) were a better fit to the solar data. We therefore conclude that the 
non-LTE corrections in AGSAK04 are overestimated, which has a significant effect on the 
concordance of the different oxygen indicators. 

AGSAK04 apphed larger downward reductions to the HM model than to the 3D and 
MARCS models, which we do not beheve to be justified. AUende Prieto et al. (2004) indicate 
that similar corrections are obtained for Kurucz and 3D models in their detailed study of 
the triplet as a function of limb darkening. This is particularly important because the 
discordance between the oxygen derived for the ID models from atomic and molecular lines 
was used as a primary argument for the superiority of the 3D models, and this discrepancy 
can be directly traced to the assignment of very large NLTE corrections to the HM model. 
By very similar logic, the internal dispersion in abundance for the atomic lines in the ID 
case arises from the assignment of large NLTE corrections to some of the hues; the internal 
agreement of the HM case is improved (and that of the 3D and MARCS cases degraded) with 
smaller NLTE effects. We illustrate this point in Figure 1. For the range of NLTE corrections 
that we consider reasonable (shaded band) the internal dispersion for the 3D models exceeds 
that of either ID model. This emphasizes the role of substantial NLTE corrections to the 
scientific conclusions of AGSAK04. 

In order to quantify this effect, we normalized the non-LTE corrections in Table 3 of 
AGSAK04 to an average for the triplet of three values: 0.16 dex (the best fit from Allende 
Prieto et al. (2004) and 0.11 dex (the mean between Holweger 2001 and Allende Prieto et 
al. (2004)), and 0.06 dex (Holweger 2001). We adopt the 0.11 dex level as our best case for 
reasons outlines below. The NLTE corrections for the other lines were scaled by the linear 
ratio of the average NLTE corrections for the triplet and the target values. We present 
revised abundance estimates for the permitted lines in AGSAK04 computed in this manner 
in Table 2. In Table 2, the first set of values contains the LTE results. The next three 
sets represent calculations where the NLTE corrections were normalized to obtain a mean 
triplet correction of 0.06 dex, 0.11 dex (our adopted mean), and 0.16 dex. The final set 
of results includes the NLTE corrections originally applied in AGSAK04. This procedure 
yields non-LTE [0/H] abundances of (8.69, 8.73, and 8.76) for 3D models, ID HM74, and 
ID MARCS respectively in the 0.16 dex case and (8.73, 8.76, 8.80) in the 0.11 dex case. 

Since the NLTE corrections are significant for the triplet, the uncertainty in these cor- 
rections is a major ingredient in the error budget. Even the reduced NLTE corrections of 
Allende Prieto et al. (2004) for the triplet are substantially larger than the corrections used 
by Holweger (2001), who found an average NLTE correction of -0.06 dex, 0.10 dex lower than 
the value reported by AGSAK04. Surprisingly, none of the authors involved commented on 
the origin of the difference. Holweger (2001) obtained average LTE and NLTE triplet abun- 
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dances of 8.78 and 8.72 respectively; his NLTE abundance is close to that obtained for the 
two ID models in AGSAK04. From AUende Prieto et al. (2004), the case with no hydrogen 
collisions was ruled out at the 3a level, and LTE models were ruled out with high confi- 
dence. However, the authors did not consider whether even lower NLTE corrections than 
their Sh = 1 case would have provided improved fits to the data. In cases such as this, 
we see no justification for simply adopting one NLTE correction (0.16 dex) over another 
pubhshed value (0.06 dex), and adopt the average of the two (0.11 dex). We note that our 
central values are close to what we would infer if we simply took the triplet alone as an 
oxygen abundance indicator and assigned a 0.16 dex NLTE correction. 

We used the dispersion in the abundances derived from individual lines as a base random 
error. A change of 0.05 dex in the triplet NLTE correction yields an average change in [0/H] 
of 0.035 dex, which we include as an additional systematic error. Finally, the log (gf) 
values from AGSAK04 are lower than previously published values by Biemont et al. (1991); 
we therefore add another 0.025 dex systematic error, following a similar error analysis by 
Mclcndcz (2004). The net effect is a total error estimate (summarized in Table 1) of 0.064 
to 0.074 dex. 



2.2.3. Infrared Molecular Lines 

The IR molecular oxygen lines have been the primary abundance indicator used in 
previous compilations of solar abundances (e.g. Grevesse & Sauval 1998), and ID model 
atmospheres yield relatively high absolute oxygen abundances. The absolute abundances 
are a strong function of the thermal structure of the model atmospheres, and the different 
thermal structure of the 3D model of AGSAK04 yields much lower predicted abundances 
than the ID models. We note that Holweger (2001) discounted CNO abundances derived 
from molecular lines because of their high temperature sensitivity. 

AGSAK04 considered two sets of OH lines: (v,r) and (r,r). The abundances predicted 
as a function of excitation potential from the 3D hydro simulations are compared with those 
from the ID Holweger- Mueller and MARCS codes in Figures 2 and 3. On these figures we 
have also indicated the average atomic abundances for each class of models. We note the 
presence of striking trends with excitation potential in the 3D and MARCS models for the 
(r,r) hues. The correspondence between the MARCS and 3D trends indicates that the origin 
of these features is common to both models, which indicates that it is a feature of the base 
atmospheric treatment rather than being induced by the convection simulation. 

In the context of ID models, trends such as those seen in the 3D models would be inter- 
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preted as a problem with the thermal structure or the assumed microturbulence. AGSAK04 
noted this trend, and claimed that it could be removed by invoking an outer atmosphere 
structure even cooler than the one predicted by the simulations. Figure 1 suggests an al- 
ternate explanation, namely that the thermal structure in the outer layers is closer to the 
hotter semi-empirical HM74 model. We will return to this point when we consider more 
recent work on CO abundances in the outer solar atmosphere. 

AGSAK04 discarded the (r,r) data for the weaker and stronger lines, in effect deriving an 
abundance from the valley in Figures 2 and 3. There is no better justification for discarding 
the high than the low points in this figure; such a procedure is not required for the fD 
models. We therefore derived average abundances using all of the features for all three 
models and both indicators; the standard deviation about the mean is an indicator of the 
quality of the fit for the individual bands. Our results are summarized in Table 1. The 
average molecular abundances were obtained with a weighted mean, but a simple averaging 
of the errors underestimates the dispersion. We therefore computed a for each band around 
the weighted molecular mean in Table 1 and averaged the (v,r) and (r,r) values to obtain 
the total error in the molecular abundances presented there. 

Melendez (2004) considered a third molecular oxygen abundance indicator, and derived 
relative abundance patterns comparable to what AGSAK04 found. We do not include this 
in Table 1 because it is not clear that systematic errors between model atmospheres codes 
can be properly accounted for in a differential analysis. If we had included the Melendez 
(2004) the mean molecular abundances would have been minimally altered for 3D and HM. 
The average would be reduced for MARCS and the internal error in the MARCS [0/H] 
would be dramatically increased. This provides further evidence that there is an underlying 
issue in the thermal structure of the MARCS model. Melendez also computed abundances 
with the same indicators as AGSAK04 for a Kurucz model atmosphere, and derived similar 
abundances as would be found for the HM74 model. Ayres et al. (2006) present evidence for 
a high solar oxygen derived from CO studies; we postpone a discussion of this interesting 
result to our conclusion, in the context of tests of the solar thermal structure. 



2.2. Jf^. Oxygen Abundance and Error Analysis 

Our overall result from the reanalysis of the AGSAK04 oxygen indicators is that the 
abundances derived from the atomic indicators are systematically increased for all models. In 
the original paper, the 3D abundance estimators were found to yield consistent abundances, 
while the ID abundances from different methods were highly discordant. This conclusion no 
longer holds when the reduced NLTE corrections inferred from limb darkening studies are 
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employed. In fact, one would obtain very similar conclusions to those presented in Table 
1 from the triplet abundance of 8.72 presented in AUende Prieto et al. (2004) for the 3D 
model (e.g. the 0.16 dex case presented in Table 2). 

Rather than simply adopting one model or another as correct, we interpret the difference 
between the HM74 and 3D abundances as evidence that the thermal structure of the Sun is 
intermediate between the two. The difference between the atomic and molecular abundances 
is roughly equal in magnitude and opposite in sign between these models; the MARCS 
model yields an intermediate abundance where the two classes of indicators give the same 
abundance, but with a larger error. The error in (Atomic-Molecular) is substantial; these 
differences are formally significant only at slightly more than la. We therefore argue that 
the mean of the derived atomic abundances (8.75) is a reasonable estimator of what one 
would obtain from a model with a thermal structure capable of reproducing the atomic and 
molecular data; one would obtain 8.74 from comparing HM and 3D, and 8.76 from MARCS 
alone (which is already internally consistent). We have a random error of 0.05 dex for the 
mean atomic abundance, but this is insufficient for a total error because of the presence 
of strong systematic differences. Adopting the consistency between atomic and molecular 
indicators as a measure of goodness of fit, la deviations could make cither the 3D model 
([0/H]=8.68) or the HM model (8.80) consistent. We treat this as a Icr = 0.06 systematic 
error, and note that it is comparable to the zero-point shift that we obtain for the atomic 
indicators relative to AGSAK04 estimated below. 

We can also examine systematic errors by comparing the AGSAK04 values with abun- 
dance estimates by other authors. These are most easily analyzed by comparing LTE abun- 
dances for the triplet. Holweger (2001) derived an average LTE triplet abundance of 8.78 for 
his standard model and 8.85 for the alternate VAL model. Biemont et al. (1991) reported 
8.84 for the triplet for a HM model and 8.78 for a MACKKL atmosphere. These should be 
compared with 8.89, 8.87, and 8.93 for the three LTE cases in AGSAK04. We note that 
the ID cases in AGSAK04 used the equivalent widths obtained with the line broadening of 
the 3D hydro models, so there will be differences between their results and those obtained 
with other ID codes. The average LTE abundance is 8.85, with a = 0.055dex. Systematic 
differences at the 0.06 dex level are thus a reasonable estimate of the current state of the art 
for oxygen abundances when estimated with different techniques. 

Adding systematic (0.06) and random (0.05) errors in quadrature, we obtain 8.75 +/- 
0.08 as our final oxygen abundance estimate for the Sun. This is less than 1.6a below 
the helioseismic abundance, and therefore we conclude that the existence of a solar oxygen 
problem has not been demonstrated with high statistical significance. 
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2.3. Cairbon Abundance Indicators 



The overall story for carbon follows a similar path to the changes in the inferred oxygen 
abundance, and the comprehensive reanalysis of AGSAB05 for carbon has a similar logical 
structure to the 2004 oxygen paper. Although both the carbon and oxygen are reduced, the 
C/0 ratio is preserved. A cool outer solar atmosphere in the 3D models yields substantially 
reduced abundances from molecular indicators, while blending features and non-LTE effects 
reduce the carbon abundance inferred from atomic features. What distinguishes the carbon 
from the oxygen case is that the NLTE effects are smaller, and as a result one might anticipate 
a smaller offset in the atomic line abundances than the change in atomic oxygen abundance 
indicators. However, this is not the actual published result; if anything, the derived carbon 
from atomic features is lower than the molecular value for all of the models presented in 
AGSAB05. Furthermore, we will demonstrate that this effect cannot be explained by any 
of the effects used to explain the differences in the comparison with prior work given by 
AGSAB05. Until the origin of this difference is understood, we therefore have to consider 
two different systematic sets of abundances for atomic features, and the best choice of model 
hinges on which set is correct. In this section, we discuss the three classes of indicators 
(forbidden and allowed atomic lines, and molecular) in turn, and as for the oxygen abundance 
synthesize our final best estimate and error in the fourth subsection. Our overall estimates 
are presented in Table 3. The top set of values represents the original AGSAB05 values for 
the different indicators. The middle set is what we would obtain with the low AGSAB05 
normalization of the atomic abundances, while the bottom set is what we obtain with the 
higher Biemont/Holweger normalization. 



2.3.1. Forbidden Lines 

The [CI] line at 8727 A was the subject of a detailed analysis by Allende Prieto, Lambert 
& Asplund (2002). They incorporated blending from a nearby Si feature to reduce the 
equivalent width attributable to carbon from 6.5 to 5.3 mA , with a corresponding reduction 
in the inferred abundance. They also employ a lower oscillator strength than previous studies. 
Unlike the case of oxygen, the carbon abundance derived from the forbidden line is almost as 
temperature sensitive as that derived from molecular features, so the atomic versus molecular 
diagnostic is less powerful for carbon than for oxygen. We adopt the AGSAB05 central 
values for our base case, but note that there may be explained systematics in the atomic 
carbon abundances in AGSAB05 which we discuss below. We include their error estimate 
for uncertainties in the equation of state (0.02 dex), but assign a larger uncertainty to the 
atomic physics (0.04 dex) in accord with the quoted theoretical uncertainties. Our principle 
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reservation on the error budget is the uncertainty in the continuum level and the contribution 
to the equivalent width of the blend from the wing of the Si feature. Their reduced permits 
only small deviations (of order 0.01 dex) in the derived carbon abundance, but the base model 
relies upon the assumption that the underlying velocity field is exact. Although the overall 
agreement with Fe (Asplund et al. 2000a) and Si (Asplund 2000) line profiles is good, it is 
not errorless. We cannot evaluate this ingredient directly, but an estimate based upon the 
mean deviation observed in clean lines would seem to be a worthwhile exercise. For the 
present, we therefore assign the same blending uncertainty of 0.03 dex adopted by AUende 
Prieto, Lambert & Asplund (2002) to obtain a total uncertainty of 0.054 dex. 



2.3.2. Permitted Atomic Lines 

AGSAB05 considered a subset of the permitted atomic features used in previous solar 
abundance studies (e.g. Biemont et al. 1993; Stuerenburg & Holweger 1990). Stuerenburg 
& Holweger (1990) found that small NLTE corrections are required for CI, and that the 
strength of the correction depends on equivalent width. They found an average of -0.05 dex; 
if restricted to the weaker lines included in AGSAB05, their average NLTE correction would 
be -0.02 dex. AGSAB05 computed Non-LTE corrections for ID models, and the MARCS 
corrections were applied to the 3D models. Hydrogen collisions were not included in the 
NLTE corrections; this resulted in larger downward abundance revisions (an average of - 
0.08 to -0.09) than Stuerenburg & Holweger (1990). AGSAB05 note that the case of carbon 
should be an analog of oxygen, and we concur. As a result, we contend that the case with hy- 
drogen collisions should be included in the base model. Both sources indicate that including 
hydrogen coUisions roughly halves the expected NLTE correction. We therefore considered 
two cases for NLTE corrections: a maximum of half the AGSAB05 value (corresponding 
to their hydrogen collision case, average -0.04 dex) and a minimum of one quarter of the 
AGSAB05 value (corresponding to the SH90 case, average of -0.02 dex). Our best value is 
the average between the two (a mean of -0.03 dex), and the error induced by uncertainties 
in NLTE corrections is 0.01 dex; adopting the AGSAB05 hydrogen collision case would only 
have changed our mean value by 0.01 dex. We applied these proportional NLTE corrections 
to the AGSAB05 LTE results for their three classes of models (middle values. Table 3). 
There is a small reduction in the dispersion (and mean trend with equivalent width) for the 
ID models and a corresponding increased for both in the 3D model; none of these features, 
however, are drastic. The internal dispersion in the permitted atomic abundances is of order 
0.03 dex. 

A more substantial issue emerges when we compare the AGSAB05 abundances with prior 
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work, and this is true even for the LTE estimates. The mean LTE abundance for the Biemont 
et al. (1993) sample is 8.56 for the hnes in common with AGSAB05; this should be compared 
with a HM LTE value of 8.48 for the latter compilation. This offset of 0.08 dex is comparable 
to the average difference between atomic and molecular abundance indicators. The mean 
difference in equivalent width and oscillator strength for the lines in common is negligible, 
and would yield an offset of less than 0.01 dex if applied under the assumption that all of 
the hnes are on the hnear part of the curve of growth. We illustrate the differences in Figure 
4, defined in the sense (Biemont- AGSAB05). In this figure we have corrected the Biemont 
abundances to the AG SAB 05 equivalent width and oscillator strengths. The differences are 
significant even for weak lines, suggesting that differences in the classical line broadening 
are probably not responsible. A similar, but smaller, effect is present in the forbidden line. 
Allende Prieto, Lambert & Asplund (2002) inferred a HM abundance of 8.48, which would 
also be obtained from Stuerenburg & Holweger (1990) when a blending correction is made to 
the equivalent width. AGSAB05 could not trace a comparable difference (0.06 dex) relative 
to the earlier work of Lambert. The only obvious source that we can derive is a note by 
Sturenburg & Holweger that they corrected their atomic abundances for the fraction of C 
tied up in CO, which could be of the right order to explain the differences. Until the origin 
of this discrepancy (which is not present for oxygen) is explained, we have to treat this as a 
systematic uncertainty in the atomic abundance scale. Abundances derived under this scale 
are the last set of values in Table 3. 



2.3.3. Molecular Lines 



We consider the same four molecular indicators that were included in AG SAB 05. They 
chose to disregard one of them (CH electronic hnes) in their derived mean abundances, on 
the grounds that they are located in a crowded portion of the spectrum and sensitive to the 
treatment of line broadening. However, the formal errors in the CH electronic abundances 
arc similar to those for the other molecular species, and as such we sec no obvious reason 
to exclude them. We do treat the CH (v,r) abundances as being more reliable, as they arc 
based on many more lines than the other diagnostics. We therefore assigned double weight 
to the CH values and single weight to both the C2 electronic and CH electronic values. As 
for oxygen, the mean was derived by a weighted average of the carbon obtained with different 
molecular indicators, and the scatter of the individual line measurements for all diagnostics 
around the adopted mean was taken as a measure of the random error. We did not include 
carbon (or oxygen) abundances derived from CO line studies, because there are complex 
correlated errors. Had we included them, the net effect would have been to increase the 
molecular abundances relative to the atomic values. 
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2.3.4- Carbon Abundance and Error Analysis 

Our final inference concerning carbon depends on wliicli atomic abundance scale is 
adopted. If we take the low scale of AGSAB05, the 3D model atomic and molecular abun- 
dances (8.40, 8.42) are closer than those for HM74 (8.45, 8.55); ID MARCS abundances are 
also consistent (8.40, 8.44). We would estimate a mean value of 8.41-8.42, with a random 
error of 0.04 dex. The HM74 average of 8.50 would be an effective 2a internal inconsistency, 
implying a 0.04 dex systematic uncertainty for a total abundance of 8.41 +/- 0.06. Adopting 
the higher scale would give pairwise results of (8.44, 8.42), (8.50, 8.55), (8.44, 8.44); all three 
models are internally consistent within the errors, and a mean abundance would be 8.47 with 
a total error of 0.05 (0.04 random, 0.03 systematic). We adopt the mean of these approaches 
(8.44), and estimate an error of 0.04 (random) and 0.04 (systematic) for a total of 0.06 dex 
when combined in quadrature. 

2.4. Nitrogen Abundance Indicators 

Our discussion of nitrogen is necessarily briefer than that of oxygen and carbon, largely 
because the published results are preliminary and incomplete. Holweger (2001) derived 
a non-LTE [N/H] = 8.0 +/- 0.11, comparable to results in previous compilations of solar 
abundances from Grevesse & Sauval (1998). The compilation of models in Asplund, Grevesse 
& Sauval (2005) yields atomic and molecular nitrogen abundance estimates of (7.85 -I-/-0.08, 
7.73 +/- 0.05), (7.97 +/- 0.08, 7.95 +/- 0.05), (7.94 +/- 0.08, 7.82 +/- 0.05) for 3D, HM, 
and MARCS respectively. The same correspondence between 3D and MARCS that was seen 
in oxygen is replicated in nitrogen, but the internal consistency in the HM model is higher 
than that in the other models. The formal significance of the disagreement in the 3D models 
is under 2a, however, so we cannot exclude the possibility that they may be consistent. We 
therefore adopt the HM result as the central value (7.96 +/- 0.06), and treat the difference 
with the 3D result (7.78 +/- 0.06) as a 2a systematic error. This yields a total uncertainty 
in [N/H] of 0.10 dex dominated by systematic uncertainties. 

2.5. The Solar Beryllium Abundance 

There is an interesting linkage between the solar O and Be abundances. In stellar 
interiors Be is destroyed at modest temperatures (of order 3.5 million K). It can therefore 
be used as a diagnostic of mixing in stars, especially in conjunction with the more fragile 
light element Li (Pinsonneault 1997). Traditional model atmospheres studies (Chmielewski 
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et al. 1975) yield a solar photospheric beryllium abundance roughly half of the meteoritic 
abundance. However, the only accessible Be feature is located in a crowded portion of the 
spectrum in the near UV, and the continuum opacity is uncertain in this regime (largely from 
the contribution of numerous weak iron lines). Since the strength of a line is a function of the 
ratio of the line to the continuous opacity, a higher photospheric Be could be derived if the 
continuous opacity background was higher than that of the model. In an important paper, 
Balachandran & Bell (1998) pointed out that nearby OH lines could be used to test the 
continuous opacity close to beryllium. They derived a UV OH lines that were too strong if 
they used the absolute oxygen abundance obtained from the IR OH lines, and interpreted this 
as evidence that the continuous opacity is underestimated in the spectral window relevant 
for Be. Similar conclusions for 3D models were obtained by Asplund (2004). Following 
Lodders (2003), we note that the uncertainties in the ad hoc corrections are substantial. 
Asplund (2004) quotes photospheric and meteoritic abundance errors of 0.09 and 0.08 dex 
respectively, implying that his zero net photospheric depletion has a la uncertainty of 0.12 
dex. Even if the Balachandran and Bell argument is entirely correct, the data sets a 2a limit 
of 0.24 dex on beryUium depletion and does not require that it be zero. There is also the 
possibihty of substantial NLTE corrections to the UV OH lines, which would reduce or even 
eliminate the requirement for a mechanism to reduce the strength of the lines. We also note 
that the value of the oxygen used by Balachandran & Bell (1998) for the IR lines (8.91) in 
the HM74 model is larger than the value derived from other molecular and atomic indicators, 
and even slightly larger than the value from the (v,r) transitions in the same model from 
the work of Asplund and collaborators. We contend that this promising approach still has 
substantial errors, including large uncertainties in the absolute solar oxygen abundance. We 
therefore beheve that the approach of Lodders (2003) is the best current picture of the degree 
of beryllium depletion in the Sun: namely, there is a substantial uncertainty in the degree 
of solar beryllium depletion, and that further work is required before powerful observational 
bounds can be used to constrain interiors calculations. 



3. Conclusions and Future Tests 

Our basic conclusion is simple: the difference between the solar CNO abundances as 
derived from model atmospheres and model interiors considerations is not statistically sig- 
nificant. The systematic errors in photospheric abundance indicators will have to be reduced 
before a "solar abundance problem" can be established (or ruled out) with confidence. How- 
ever, the disagreement between the solar thermal structure and that of the simulations would 
favor the higher abundance scale, and there is some recently published evidence to that ef- 
fect. If this is confirmed, it switches the nature of the problem from being a question of 
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the correct abundance scale to a question of the uncertainties in numerical convection sim- 
ulations. We begin with a synthesis and explanation of our findings. We then divide our 
conclusions into two parts. We recommend steps to more firmly establish the photospheric 
abundance scale, and contend that accurate solar abundances require tests of the thermal 
structure of the models and the magnitude of non-LTE abundance corrections. In our final 
subsection we then gather together evidence that the atmospheric abundance scale problem 
may be tied to the limited resolution in the convection simulations or errors in the underlying 
model atmosphere treatment. The consequences for the solar beryllium abundance, which 
is a useful diagnostic of internal mixing, are also explored. 

The two main justifications for the superiority of the 3D hydro atmospheres are the 
treatment of line broadening and the inclusion of granulation. Both of these represent gen- 
uine improvements in the atmospheric physics. However, neither of these effects is actually 
primarily responsible for the difference in the solar abundance scale. Many of the abun- 
dance indicators are insensitive to the effective microturbulence. If temperature fluctuations 
are imposed on a semi-empirical Holweger-Mueller atmosphere, the resulting granulation 
corrections are usually smaller than the 3D convection effects reported by Asplund and col- 
laborators, and frequently opposite in sign (Holwcgcr 2001). The main driver behind the 
systematic reductions in abundance derived from the 3D models is a theoretically predicted 
change in the thermal structure, coupled with large assumed non-LTE corrections for atomic 
features. Neither of these changes is directly supported by observational tests. Instead, the 
argument for the superiority of the abundances derived from the newer model atmospheres is 
an indirect one, focused on the concordance of abundances derived from different indicators. 

A consistent chain of logic emerges from the comprehensive studies of oxygen (AGSAK04) 
and carbon (AGSAB05). Classical LTE model atmospheres tend to yield internally consis- 
tent, and high, carbon and oxygen abundances for atomic and molecular indicators. The 
application of a different thermal structure in the 3D hydro atmospheres drastically reduces 
the abundances inferred from highly temperature sensitive molecular indicators, but has a 
smaller effect on atomic features. Large NLTE corrections are then applied to the abun- 
dances derived from permitted atomic features for both ID and 3D models. The net result is 
that the abundance estimates from ID models become internally inconsistent (atomic indi- 
cators yield lower abundances than molecular ones) , while abundances derived from the 3D 
models are internally consistent. The abundances derived from forbidden lines are insensitive 
to NLTE effects, but they are reduced in the newer generation of models by the inclusion of 
blending features. As a secondary argument, the fits to individual indicators are argued to 
be superior in the 3D models when compared to the fits to individual indicators in the ID 
models. This approach is appealing on the surface, but when examined in detail the picture 
is decidedly more ambiguous. If anything, the hints from the data would lean towards the 
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opposite conclusion. 

The abundances derived from forbidden lines have the smallest systematic errors, but 
errors in both the theoretical oscillator strengths and the treatment of blending features result 
in non-negligible random errors. More to the point, the internal consistency of abundances 
derived from forbidden and molecular lines is actually similar in the 3D and ID cases. Prom 
Table 1, the forbidden and molecular oxygen abundances are (8.71, 8.64) for 3D and (8.77, 
8.84) for the HM; the differences are identical. Given the errors, neither discrepancy is 
statistically significant with high confidence. 

The abundances reported for permitted atomic features in AGSAK04 and AGSAB05 
are significantly lower for ID models than the corresponding molecular abundances, while 
the reported 3D results are in agreement. In the case of oxygen, this rests completely on the 
assignment of large NLTE corrections. These corrections were obtained under the assumption 
that hydrogen collisions were unimportant. Detailed studies of the response of the triplet 
to limb-darkening indicate that models including hydrogen collisions are favored, and the 
inferred NLTE corrections decrease. As a result, the internal consistency of the oxygen 
indicators is comparable for the different classes of atmospheres. Nitrogen is consistent for 
HM74 models and inconsistent (but at less than 2a) for the 3D case. In the case of carbon, 
the situation is made more complex by significant zero-point offsets between earlier studies of 
carbon abundances that arc not explained. Again, the assignment of larger NLTE corrections 
is uncertain (and, unlike the case of oxygen, not directly tested against limb-darkening data). 
A clean distinction between models on the basis of consistency is not obtained. However, 
the 3D models do yield different molecular and atomic abundances for both N and O, and 
might also do so for C. 

One might then hope to find distinct differences in the quality of the fits to different 
molecular indicators. The usual patterns, unfortunately, manifest themselves as simple zero- 
point shifts. For every case where there are issues with the ID models (e.g. small trends 
with excitation potential in the [0/H] derived from (v,r) OH transitions in the HM model) 
there are comparable or larger effects for the 3D models (e.g. substantial trends in the [0/H] 
derived from (r,r) OH transitions). In a recent preprint, Scott et al. (2006) examined CO 
indicators, and the resulting pattern is illustrative. The 3D models yielded similar results 
for two of the three features studied, while the ID models performed better in a different 
pair of indicators. The C^^/C^^ ratio from the ID models ranges from 69 to 84, while the 
same ratio for the 3D models ranges from 83 to 108. These values should be contrasted with 
the expected terrestrial ratio of 89. Scott et al. (2006) choose comparisons that favor the 
3D models, while an advocate of the traditional models might reasonably stress the other 
cases. In our view, the best choice of models is not clearly distinguishable from the CNO 
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abundance studies. We recommend caution when extrapolating these model results to other 
stars, where the differential effects can be even more drastic. 

3.1. Establishing the Absolute Photospheric Abundance Scale 

The single most important test that is required for atmospheres theory is a discriminant 
between the different proposed thermal structures of the solar atmosphere. The recent paper 
by Ayres et al. (2006) makes an important contribution by making direct comparisons of 
solar data with the thermal properties of the simulations. They present evidence that the 
solar center-to-limb variations in continuum flux arc inconsistent with the predictions of 
the 3D hydro simulations. They also note that the predicted magnitude of fluctuations in 
the upper atmosphere from the simulations appears to be larger than the observed pattern. 
Ayres et al. then construct an empirical model of the atmosphere and derive a high oxygen 
abundance (8.85) from CO molecular features under the assumption of a fixed C/0 ratio. 
In retrospect this conclusion is not surprising. The HM model is not a purely theoretical 
exercise; it was constructed to reproduce the mapping of the source function as a function of 
optical depth inferred from limb darkening studies of continuum flux and strong lines (see also 
Allcndc Pricto et al. 1998). The relative trends we have inferred from atomic and molecular 
abundance indicators support the conclusions of Ayres et al., but the current errors make 
our evidence in this matter suggestive but not conclusive. It would also be highly beneficial 
to repeat the HM exercise with the full 3D models as opposed to the restricted form of them 
that Ayres et al. (2006) had available to them. 

Ultimately, the absolute accuracy of photospheric abundances is directly tied to the 
absolute accuracy of the thermal structure. This suggests that an approach similar to that 
of Steffen & Holweger (2002) may be the optimal one. In their paper they examined the 
impact of temperature fluctuations around an assumed mean empirical thermal structure, 
which in their case was the HM74 model. Interestingly, the abundance corrections that they 
derive would act in the sense of increasing the concordance between abundance indicators. 
Oxygen abundances from atomic indicators would be slightly increased; although they did 
not consider molecular features directly, the net effect would certainly have the same sign as 
that obtained from 3D hydro models, namely a decrease in the inferred abundance. In such 
a differential approach, deviations between the mean structure of the simulations and the 
empirical data would be used as guidance concerning the underlying physics. In contrast, the 
3D model abundances assume that the ab initio profile is correct. A similar approach could 
be employed for the velocity field that replaces the microturbulence and macroturbulence in 
traditional ID atmospheres. 
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A second ingredient that must be tested empirically, rather than by theoretical assertion, 
is the magnitude of NLTE corrections. The available evidence suggests that NLTE correc- 
tions are in general small for the Sun, but for the level of precision required in the absolute 
abundance scale these small corrections are significant. Studies of different spectral features 
yield different conclusions about the physical model employed in NLTE studies. This im- 
plies that there are significant uncertainties in absolute theoretical calculations. Fortunately, 
NLTE corrections can be constrained by the response of line strength to limb darkening in the 
Sun. It should be possible to develop improved theoretical models with a sufficient database 
of information developed in this fashion. One other stringent test of NLTE effects may be to 
focus on the species whose relative abundances can be reliably inferred from meteoritic data. 
For example, NLTE effects may be significant for iron (Shchukina & Trujillo Bueno 2001) 
but less so for Si (Wcdemeyer 2001). Holweger (2001) noted that there may be a conflict 
between the photosphcric and meteoritic Fe/Si ratio, albeit one of marginal significance. A 
similar situation may exist for Na (Asplund, Grevesse & Sauval 2005). 

Another tractable problem is the absolute error for the forbidden C and O lines. In 
these cases, uncertainties in the line profiles and continuum levels should be included. Better 
atomic data (such as oscillator strengths for both the lines and the blending features) would 
also be useful. 

The accuracy of the theoretically predicted turbulent velocity field as a function of 
optical depth should also be subjected to a more rigorous analysis. Scott et al. (2006) 
present evidence that the generation of simulations used for the abundance analysis yielded 
poor fits to the line bisectors of CO lines. Higher resolution simulations gave better line 
profile fits, but for (unspecified) unrealistically high C/0 abundances. The higher resolution 
simulations were not employed in the CO abundance analysis in that paper. It is worth 
keeping in mind that line profiles are integral quantities, and as a result the uniqueness of 
the solutions is not established by individual cases of good fits. This is particularly true when 
the abundance itself is treated as a free parameter. It would be extremely useful if future 
papers on abundances derived using numerical simulations illustrated individual line fits, as 
well as quantifying the actual impact of the "effective microturbulence" on the abundance 
estimates. 

It is useful to separate out the impact of velocity broadening from the effect of gran- 
ulation and temperature gradient changes. This can be done by using the mean thermal 
structure and temperature fluctuations from the simulations and a more traditional mi- 
cro/macroturbulence model to infer abundances, and comparing the results with the full 3D 
models. Scott et al. (2006) constructed such a test case (their IDAV model), and found only 
small abundance offsets, of order 0.04 dex for oxygen derived from IR OH fines. They also 
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inferred carbon abundances from CO; in this case O was held fixed and the carbon was ad- 
justed to fit different molecular indicators. The deviations in the derived carbon abundances 
relative to the 3D case ranged from small (0.01 dcx for the LE lines) to modest (0.06 dex for 
the weak Az/ = 1 lines) to large (0.14 dex for the Az/ = 2 lines). These deviations may ex- 
plain the changes in excitation potential that Ayres et al. (2006) needed to obtain consistent 
abundances within a ID framework. This exercise implies that the impact of the improved 
microphysics varies substantially for different indicators, and is worth quantifying across 
the board. An alternate exercise (using the revised velocity field and relative temperature 
fiuctuations while adopting a HM74 mean thermal structure) might also be illuminating. 



3.2. Uncertainties in Numerical Convection Simulations 

First-principles theoretical model atmosphere calculations have undeniable strengths. 
The ability to naturally reproduce line widths and include granulation is a powerful addition 
to our ability to reliably interpret stellar and solar spectra. The principal difficulty with such 
models is that errors in the input physics generate absolute errors in the inferred atmospheric 
structure that cannot be calibrated away in the absence of explicit free parameters. This 
phenomenon is the major reason why numerical convection simulations have not replaced 
the simple mixing length theory in stellar interiors calculations. Interiors models that can 
reproduce observed stellar radii are simply more useful for most purposes than models with 
a better physical treatment of convection that fail to do so. 

Before the results from such models are adopted as the new abundance standard, it will 
be necessary to perform an extensive theoretical error analysis and to compare the models 
with the strongest observational constraints. We believe that accurate solar abundance 
calculations must reproduce the observed solar thermal structure, and from the Ayres et 
al. (2006) paper the Asplund models employing numerical convection simulations appear to 
yield a temperature gradient steeper than the real Sun. This could be caused by errors in the 
background (ID) stellar atmospheres treatment; for example, uncertainties in the equation 
of state and continuous opacities will induce absolute errors in the thermal structure. An 
approach similar to that employed in interiors models would be useful for assessing the 
uncertainties in the thermal structure and abundance predictions, and this should be included 
in the error budget for abundances. 

It is more likely, however, that the major error source in 3D hydro model atmospheres 
is related to uncertainties in the numerical convection simulations. The approximations in 
hydro simulations of giant planet atmospheres are have been demonstrated to be strongly 
affected by the quality of the assumed physics (Evonak & Glatzmaier 2004). Zhang & 
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Schubert (2006) also provides a good summary of the uncertainties in the related prob- 
lem of terrestrial and solar dynamo models. Another phenomenon that could be related 
is the issue of convective overshooting below surface convection zones. Numerical simula- 
tions have tended to favor extensive overshooting, and the early models had a substantial 
nearly adiabatic overshoot region, in conflict with the stringent limits set by seismology (less 
than O.OSifp). More recent 3D Brummell et al. (2002) and 2D (Rogers & Glatzmaier 2005, 
2006) calculations found that the fiUing factor for plumes is smaller than previously thought, 
which led to an overestimate in earlier models of the changes induced by overshooting in 
the thermal structure. The newer simulations predict strongly subadiabatic overshooting 
(effectively, overmixing without changing the thermal structure), which is consistent with 
the seismic limits. However, they still produce a substantial mixed region below the surface 
convection zone of order OAHp. Since even a small overmixing of 0.05Hp drastically increases 
pre-MS lithium depletion (Pinsonneault 1997), which is already too efficient relative to stel- 
lar data (Piau & Turck-Chiczc 2002), it is likely that even this reduced overshooting is too 
large to be compatible with stellar constraints. We argue that there is a common pattern 
in both "undershooting" and "overshooting" above and below convective regions. In both 
cases, the numerical simulations may be overestimating the degree of mixing and the impact 
on the thermal structure of convection outside the formal bounds set by the Schwartzschild 
criterion. 

There are two plausible error sources that should be investigated. The treatment of heat 
transfer in the atmosphere convection simulations is necessarily simplified, and this may be 
leading to an artificial inhibition in energy transport between turbulent cells projected into 
the radiative atmosphere and their surroundings. Resolution effects, however, may be even 
more important. Even the highest resolution simulations available today are many orders of 
magnitude away from being able to reproduce the characteristic Reynolds numbers in the 
Sun. Scott ct al. (2006) found significant changes in line bisectors for the outer layers of 
their solar model when they increased their resolution, and these changes were in the sense 
of reducing the temperature contrast in the upper atmosphere and improving the shape of 
the bisectors relative to data. Numerical tests with substantially increased resolution may 
shed some interesting light on the sensitivity of the predictions to the underlying numerics; 
2D convection simulations may be useful in this regard. We are optimistic that the net effect 
of such testing will be a greatly improved understanding of the strengths and weaknesses 
of theoretical atmospheres models, just as we are confident that the net result of the solar 
abundance controversy will be a far more secure knowledge of stellar abundances. 

We would like to thank Martin Asplund for providing tables of the abundances derived 
from molecular oxygen abundance indicators. We would also like to thank Don Terndrup, 
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Indicator 


Mean 


Error 


Mean 


Error 


Mean 


error 




OLJ 




"RMnzL 








[0 I] 


8.680 


0.010 


8.760 


0.010 


8.720 


0.010 


I 


8.639 


0.021 


8.644 


0.065 


8.716 


0.031 


OH (v,r) 


8.610 


0.030 


8.870 


0.030 


8.740 


0.030 


OH(r,r) 


8.650 


0.020 


8.820 


0.010 


8.830 


0.030 


Atomic 


8.660 


n/a 


8.702 


n/a 


8.718 


n/a 


Moleculcir 


8.630 


n/a 


8.845 


n/a 


8.785 


n/a 


Atom-Mol 


0.030 


n/a 


-0.143 


n/a 


-0.067 


n/a 




Rev3D 




RevHM 




RevMARCS 




[0 I] 


8.700 


0.055 


8.780 


0.055 


8.740 


0.055 


I 


8.721 


0.074 


8.742 


0.065 


8.794 


0.064 


OH (v,r) 


8.612 


0.028 


8.865 


0.030 


8.740 


0.027 


OH(r,r) 


8.674 


0.037 


8.826 


0.017 


8.850 


0.050 


Atomic 


8.707 


0.044 


8.764 


0.042 


8.763 


0.042 


Moleculcir 


8.635 


0.048 


8.835 


0.031 


8.765 


0.082 


Atom-Mol 


0.073 


0.065 


-0.071 


0.052 


-0.002 


0.092 



Table 1: Comparison between the O abundances derived from several indicators. The upper 
part of the tables recals the value obtained by Asplund, Grevesse & Sauval (2005) while the 
bottom part summarize the our reanalysis (see text for details). 



-31 - 



Indicator 


NLTE triplet correction 


Wavelength (A) 


3D 


HM04 


M04 


I 


0.000 


6158.1 


8.650 


8.820 


8.800 


I 




7771.9 


8.910 


8.890 


8.950 


I 




7774.2 


8.890 


8.870 


8.940 


I 




7775.4 


8.860 


8.880 


8.910 


I 




8446.7 


8.680 


8.690 


8.750 


I 




9266.01 


8.730 


8.790 


8.770 


Mean 






8.803 


8.840 


8.872 


Dispersion 






0.105 


0.071 


0.084 


I 


0.060 


6158.1 


8.642 


8.806 


8.792 


I 




7771.9 


8.836 


8.811 


8.885 


I 




7774.2 


8.825 


8.796 


8.877 


I 




7775.4 


8.805 


8.809 


8.855 


I 




8446.7 


8.658 


8.660 


8.728 


I 




9266.01 


8.708 


8.760 


8.748 


Mean 






8.758 


8.787 


8.829 


Dispersion 






0.080 


0.055 


0.064 


I 


0.110 


6158.1 


8.635 


8.795 


8.785 


I 




7771.9 


8.775 


8.745 


8.830 


I 




7774.2 


8.770 


8.735 


8.825 


I 




7775.4 


8.760 


8.750 


8.810 


I 




8446.7 


8.640 


8.635 


8.710 


I 




9266.01 


8.690 


8.735 


8.730 


Mean 






8.721 


8.742 


8.794 


Dispersion 






0.060 


0.049 


0.048 


I 


0.160 


6158.1 


8.628 


8.784 


8.778 


I 




7771.9 


8.714 


8.679 


8.775 


I 




7774.2 


8.715 


8.674 


8.773 


I 




7775.4 


8.715 


8.691 


8.765 


I 




8446.7 


8.622 


8.610 


8.692 


I 




9266.01 


8.672 


8.710 


8.712 


Mean 






8.684 


8.697 


8.759 


Dispersion 






0.041 


0.052 


0.035 


I 


variable 


6158.1 


8.620 


8.770 


8.770 


I 




7771.9 


8.640 


8.600 


8.710 


I 




7774.2 


8.650 


8.600 


8.710 


I 




7775.4 


8.660 


8.620 


8.710 


I 




8446.7 


8.600 


8.580 


8.670 


I 




9266.01 


8.650 


8.680 


8.690 


Mean 






8.639 


8.644 


8.716 


Dispersion 






0.021 


0.065 


0.031 



Table 2: Oxygen abundance derived from the O I hnes using different NLTE corrections. 
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Indicator 




tirror 


iviean 


tirror 




error 




3D 




HM 




M 




[CI] 


8.390 


0.040 


8.450 


0.040 


8.400 


0.040 


C I 


8.360 


0.030 


8.390 


0.030 


8.350 


0.030 


CH - vib-rot 


8.380 


0.040 


8.530 


0.040 


8.420 


0.040 


C2 electronic 


8.440 


0.030 


8.530 


0.030 


8.460 


0.030 


CH electronic 


8.450 


0.040 


8.590 


0.040 


8.440 


0.040 


Average 


8.394 


0.036 


8.483 


0.036 


8.405 


0.036 


Deviation 


0.034 




0.068 




0.046 




ATM 


8.375 




8.420 




8.375 




MOL 


8.413 




8.545 




8.435 




ATM-MOL 


-0.038 




-0.125 




-0.060 




Corrected Low Atomic Scale 


3D 




HM 




M 




[CI] 


8.390 


0.040 


8.450 


0.040 


8.400 


0.040 


C I 


8.406 


0.042 


8.441 


0.035 


8.398 


0.031 


CH - vib-rot 


8.380 


0.051 


8.530 


0.039 


8.420 


0.039 


C2 electronic 


8.440 


0.037 


8.530 


0.030 


8.460 


0.033 


CH electronic 


8.450 


0.057 


8.590 


0.058 


8.440 


0.039 


Average 


8.405 


0.045 


8.495 


0.039 


8.417 


0.036 


Deviation 


0.026 




0.049 




0.029 




ATM 


8.398 




8.445 




8.399 




MOL 


8.413 




8.545 




8.435 




ATM-MOL 


-0.014 




-0.100 




-0.036 




Corrected high Atomic Scale 


3D 




HM 




M 




[CI] 


8.425 


0.040 


8.485 


0.040 


8.435 


0.040 


C I 


8.486 


0.042 


8.521 


0.035 


8.478 


0.031 


CH - vib-rot 


8.380 


0.051 


8.530 


0.039 


8.420 


0.039 


C2 electronic 


8.440 


0.037 


8.530 


0.030 


8.460 


0.033 


CH electronic 


8.450 


0.057 


8.590 


0.058 


8.440 


0.039 


Average 


8.434 


0.045 


8.524 


0.039 


8.446 


0.036 


Deviation 


0.044 




0.021 




0.026 




ATM 


8.456 




8.503 




8.457 




MOL 


8.413 




8.545 




8.435 




ATM-MOL 


0.043 




-0.042 




0.022 





Table 3: Comparison of the different C abundances derived from the indicators used by 
Asplund et al. (2005) when differents set of correction are applyed. The top panel corresponds 
to the results derived in Asplund et al. (2005), the middle panel correspond to the revised 
value when we adopt another NLTE correction that is smaller than those of Asplund et al. 
(2005) and the bottom panel includes a shifth of -1-0.035 dex and -1-0.08 dex in the forbiden 
and allowed lines respectively (see text for detail) . 
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0.1 0.2 0.3 

Triplet NLTE corrections 



Fig. 1. — The internal scatter in oxygen inferred from permitted atomic lines as a function 
of the NLTE correction applied to the triplet. The solid line is the 3D model, the dotted 
hne is the ID HM74 model, and the dash-dotted line is the ID MARCS model. AGSAK04 
adopted the values on the right side. The band denotes the limiting cases discussed in the 
text. The 3D models are less internally consistent than the ID models for NLTE corrections 
in this range. 
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Fig. 2. — The oxygen abundances derived from infrared OH lines in the ID HM74 and 3D 
AGSAK04 atmospheres are compared as a function of excitation potential The left panel 
represents (r,r) transitions, while the right panel represents (v,r) transitions. Results from 
the 3D model atmospheres are solid circles; results from the HM74 semi-empirical model 
atmosphere are open circles. The solid and dashed lines indicate the mean abundances that 
would be derived from atomic features in the 3D and HM74 models. Comparable discrepan- 
cies between molecular and atomic abundances are present for both classes of models, and 
the 3D models exhibit striking trends with excitation potential not seen in the HM74 case. 
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Fig. 3. — The oxygen abundances derived from infrared OH lines in the ID MARCS and 3D 
AGSAK04 atmospheres are compared as a function of excitation potential The left panel 
represents (r,r) transitions, while the right panel represents (v,r) transitions. Results from 
the 3D model atmospheres are solid circles; results from the MARCS theoretical ID model 
atmosphere are open circles. The solid and dashed lines indicate the mean abundances that 
would be derived from atomic features in the 3D and MARCS models. The MARCS model 
is more internally consistent than the 3D case, but similar trends with excitation potential 
are seen in both. Since the underlying model atmosphere code is similar, we conclude that 
the origin of these trends is in the basic model atmospheres rather than being induced by 
the treatment of convection. 
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Fig. 4. — Top panel: Comparison of the derived oxygen abundance from the atomic hnes 
common to Asplund, Grevesse & Sauval (2005) (filled circles) and Biemont et al. (1993) 
(empty circles) as a function of the equivalent width. The data from Biemont ct al. (1993) 
have been corrected for the difference in log (gf) and in equivalent width, bottom panel: 
Difference between these abundances as a function of equivalent width. 



